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1  Introduction  and  motivation 

It  is  well-known  that  there  can  be  a  significant  difference  in  behavior  in  flight  between 
liquid-filled  and  solid-filled  projectiles.  The  difference  is  caused  by  the  motion  of  the  liquid 
inside  the  spinning  projectile.  This  motion  causes  forces  to  act  on  the  projectile  that  can 
ultimately  cause  the  flight  of  the  projectile  to  be  prematurely  terminated  by  instability. 
The  initial  motion  of  the  projectile  necessarily  causes  the  fluid  motion  in  the  cylinder  to  be 
time-dependent;  later  it  can  be  assumed  that  the  flow  is  steady.  The  destabilizing  motion 
of  the  projectile  induced  by  the  forces  exerted  by  the  fluid  is  of  small  amplitude  in  its 
initial  stages. 

Several  methods  are  available  to  find  the  fluid  motion  in  a  cylindrical  container  per¬ 
forming  small  amplitude  oscillations;  here  we  mention  two  of  these  approaches.  First,  a 
finite  difference  approach  to  the  equations  of  motion  can  be  used1 .  In  this  formulation,  the 
Navier-Stokes  equations  are  marched  forward  in  time  until  an  equilibrium  state  is  reached. 
This  is  not  a  time-accurate  method,  however.  The  method  is  most  economical  at  relatively 
small  Reynolds  numbers  and  prohibitively  expensive2  at  Reynolds  numbers  of  the  order  of 
10,000. 

Secondly,  a  spatial  eigenvalue  procedure  has  been  developed  by  Hall,  Sedney,  and  - 
Gerber3.  This  procedure  can  be  used  to  determine  inexpensively  the  forced  fluid  mo¬ 
tion  at  Reynolds  numbers  up  to  about  2000.  The  method  produces  results  in  excellent 
agreement  with  experimental  observations  and  agrees  with  the  finite  difference  results  in 
most  situations.  A  solution  is  obtained  with  this  method  by  expressing  the  velocity  and 
fluid  pressure  fields  in  terms  of  the  eigenfunctions  that  describe  the  linear  instability  of 
solid-body  rotation.  However,  because  we  are  concerned  with  flows  in  cylinders  of  finite 
length,  we  can  allow  the  axial  wavenumber  of  the  perturbations  to  be  complex  while  the 
frequency  of  the  perturbations  is  real.  The  eigenvalue  problem  associated  with  this  type  of 
perturbation  has  an  infinite  number  of  eigenvalues,  which  can  be  numbered  in  some  way, 
for  example,  by  counting  the  number  of  internal  zeros  of  the  eigenfunctions.  Essentially  it 
is  found  that  eigenvalues  can  be  naturally  split  into  groups  of  three  with  the  eigenvalues 
within  each  group  now  being  ordered  on  the  basis  of  the  number  of  zeros  of  the  eigenfunc¬ 
tions.  The  major  computational  task  associated  with  the  spatial  eigenvalue  approach  is 


k 
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to  determine  the  eigenvalues.  Most  of  the  calculations  performed  by  the  present  authors 
were  carried  out  by  using  initial  guesses  for  the  eigenvalues  from  interpolation  in  a  table  of 
eigenvalues4  obtained  by  extrapolation  from  the  eigenvalues  at  a  Reynolds  number  of  10. 
In  order  to  perform  calculations  over  a  wider  range  of  parameter  values  it  was  necessary 
to  find  an  alternative  (and  quick)  way  of  obtaining  initial  guesses  for  almost  all  of  the 
eigenvalues.  This  report  demonstrates  the  most  useful  such  method,  which  is  based  on  the 
WKB  method. 

In  this  report  we  explain  the  implementation  of  the  WKB  method  to  find  the  eigen¬ 
value  spectrum  for  solid-body  rotation.  The  reader  is  referred  to  References  3  and  4 
for  discussions  of  the  spatial  eigenvalue  approach.  Here  the  approach  is  outlined  for  the 
completely  filled  cylinder  problem;  the  extension  of  the  approach  to  the  ‘inner-rod’  or 
‘partially-filled’  case  is  straightforward. 

2  The  spectrum  at  high  Reynolds  numbers 

As  mentioned  above,  the  motivation  for  this  work  comes  from  our  previous  study  of 
forced  motions  in  finite  rotating/nutating  cylinders3.  We  first  present  a  brief  description 
of  the  physical  problem  that  motivates  the  present  study.  Consider  a  cylinder  of  radius 
a  and  height  2c,  rotating  about  its  axis  with  angular  velocity  Cl.  The  cylinder  is  filled 
with  an  incompressible  fluid  of  density  p  and  kinematic  viscosity  v.  After  transient  effects 
(spin- up)  have  died  out,  the  fluid  motion  in  the  absence  of  gravity  is  solid  body  rotation. 
If  the  cylinder  is  made  to  cone  at  constant  rate  r  with  respect  to  an  inertial  reference 
frame,  the  motion  is  a  perturbation  on  solid  body  rotation.  In  coning  motion,  the  axis  of 
the  cylinder  moves  on  a  cone;  for  convenience,  the  vertex  of  the  cone  is  taken  to  be  at  the 
geometrical  center  of  the  cylinder.  The  angle  between  the  cylinder  and  the  cone  axis  is 
Kq.  We  define  the  Reynolds  number  Re  by 

(Cl  +  r  cosKo)a 2 

Re  —  - 

v 

and  the  nondimensional  perturbation  frequency  by 


/  = 


r 

n  +  r  cos  Kq 
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For  the  linearized  problem,  cosKo  =  1.  Time,  length,  velocity  and  pressure  are  made 
nondimensional  by  (ft  +  r)"1,  a,  (ft  +  r)a  and  pa2 (ft  +  r)2,  respectively.  Our  starting 
point,  therefore,  is  the  linearized  stability  equations  for  a  three-dimensional  time-dependent 
disturbance  to  solid  body  rotation  in  a  circular  cylinder.  We  restrict  our  attention  to  the 
mode  with  azimuthal  wavenumber  unity  and  axial  wavenumber  k.  This  wavenumber  is 
taken  to  be  complex  while  the  frequency  of  the  disturbance  is  real;  notice  that  at  finite 
values  of  the  Reynolds  number,  solid  body  rotation  is  stable  so  that  no  real  eigenvalues  k 
can  exist  .  At  infinite  Reynolds  numbers  the  flow  is  neutrally  stable  so  that  for  high  values 
of  the  Reynolds  number  we  expect  to  find  eigenvalues  k  with  small  imaginary  part.  For 
the  inertial  reference  frame,  the  appropriate  stability  equations  are: 

[i?e-1(Ai  —  r  2)  -  iM)u  +  (  2  + — — J  v  -  pr  =  0 

[Re-^Aj  -  r“2)  -  iM\v  -  ^2  +  2-~  )  u  +  ~  =  0  (2.1  a,b,c,d) 

[Re-1(Ai)  —  iM]w  +  kp  —  0 
(ru)r  —  iv  -f  krw  =  0 


Here  r  denotes  the  radial  variable  with  the  cylinder  located  at  r  =  1  while  u,t>,  u>  denote 
the  radial,  azimuthal,  and  axial  perturbation  velocities  and  p  denotes  the  perturbation 
pressure.  The  operator  Ai  is  given  by 


<P  Id  n2  1  , 

l~dr2  +  rdr  ^  +  r2  ^ 


A,  = 


(2.2) 


and  M  =  f  —  1  with  /  equal  to  the  perturbation  frequency.  If  we  are  concerned  with  the 
case  where  there  is  no  inner  cylinder  then  (2.1)  must  be  solved  for  k  =  k(Re,  f)  such  that 


u  —  iv  =  w  =  p  =  0,  r  =  0 
u  =  v  =  w  =  0,  r  =  l 


(2.3a,  b) 


The  first  of  these  conditions  ensures  that  the  velocity  and  pressure  fields  of  the  disturbance 
are  regular  at  the  origin  and  the  second  condition  is  the  no-slip  condition  for  a  viscous  fluid. 
Notice  that  if  we  were  investigating  the  ‘inner  rod’  problem  or  the  ‘two-fluid  ’problem, 
then  the  inner  boundary  conditions  would  be  altered  and  applied  instead  at  a  finite  value 
of  r. 
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Our  aim  therefore  is  to  determine  the  infinite  spectrum  of  eigenvalues  at  a  given  value 
of  /  at  high  values  of  Re.  From  our  previous  experience  (e.g.,  see  Reference  3)  we  know 
that  for  positive  values  of  the  real  part  of  k,  there  are  three  distinct  branches  of  eigenvalues. 
Two  of  these  branches,  I  and  III,  are  in  the  first  quadrant  while  a  third  branch,  II,  is  in  the 
fourth  quadrant.  We  know  that  on  any  one  of  these  branches  the  asymptotic  expansion  of 
any  eigenvalue  takes  the  form 

k  =  ko  -f  Re~ 5  k\  +  Re~l  k%  +  •  •  •  (2.4a) 

or  the  form 

k  =  Rel>k0  +  Re~h1+--.  (2.46) 

The  eigenfunction  associated  with  ko  in  (2.4a)  can  be  expressed  in  terms  of  Bessel  functions 
and  ki  defends  on  the  mode  number  m  (basically  the  number  of  zeros)  of  the  eigenfunction. 
In  fact,  for  large  m  it  can  be  shown  that  k\  ~  m,  so  that  the  first  two  terms  in  (2.4a) 
will  formally  break  down  when  m  ~  Re  ? .  In  this  situation  the  radial  derivatives  will 
formally  be  of  size  Re$  so  that  the  basic  assumptions  leading  to  (2.4a)  no  longer  apply.  It 
is  this  regime  that  we  will  now  concentrate  on  since  our  previous  experience  has  shown  that 
(2.4a,b)  are  only  useful  in  predicting  accurate  approximations  to  the  first  few  eigenvalues 
on  the  different  branches. 

For  modes  with  wavenumbers  of  size  0(Re  s )  the  disturbance  varies  on  a  length  scale 
0(/?e~  3 )  over  most  of  the  flowfield;  we  anticipate  this  structure  by  seeking  WKB  type 
solutions  of  (2.1)  with 

/r  oo 

0(f)dr]  ^(un,un,u;n,i?e'T'pn)i?e”* .  (2.5) 

n=0 

At  this  stage  the  viscous  derivatives  in  the  radial  direction  are  of  size  0 (Re)  while  diffusion 
in  the  axial  direction  will  be  of  size  0(k2).  Thus  in  order  to  balance  the  diffusion  of  vorticity 
in  the  radial  and  axial  directions  we  must  take  k  —  0(iZe*)  so  that  (2.4a, b)  must  now  be 
replaced  by  the  expansion 

k  =  VRe.  k0  -f  *i  4-  3 — k2  -f - .  (2.6) 

wRe 
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The  WKB  phase  function  9  appearing  in  (2.5)  is  determined  by  substitution  for  {u,  v,  tu,p) 
from  (2.5)  into  (2.1)  and  retaining  only  the  dominant  terms  in  the  limit  Re  ->  oo.  We 
obtain 


-\92  +kl  +  iM]u0  +  2uo  =  tfljj 


4-0 


[92  -f^o+  *-M]u0  +  2uq  =  0 
-[02  4-  kl  4-  4-  kop^  =  0 

i9uo  +  koWQ  =  0, 

and  these  equations  have  a  consistent  solution  if 


(2.66) 


(92  4-  kl  4-  iM)2(92  +  kl)  +  4kl  =  0, 


(2.7) 


so  there  are  six  possible  values  for  6  for  a  given  value  of  k.  If  the  expansion  (2.5)  is  to  be 
valid  over  a  0(  1 )  length  in  r  then  only  real  values  of  9  are  acceptable.  Thus  we  insist  that 
(2.7j  has  two  real  roots,  in  which  case  it  can  be  shown  that  the  other  four  roots  are  complex. 
Without  any  loss  of  generality  we  can  suppose  that  ±9i  are  the  two  real  roots  of  (2.7)  while 
the  remaining  (complex)  roots  are  denoted  by  ±92,±93.  Hence  our  expansion  (2.5)  must 
now  be  written  down  as  a  sum  over  the  two  acceptable  real  solutions  of  (2.7);  alternatively 
these  exponential  solutions  can  be  combined  in  terms  of  trigonometric  functions. 

At  next  order  the  linear  equations  to  determine  (ux ,  v j , )  are  obtained  but  are 
found  to  be  forced  by  the  zeroth  order  function  their  derivatives  with 

respect  to  the  slow  scale  r.  The  consistency  of  these  equations  leads  to  a  first  order 
differential  equation  for  Mo  that  can  be  solved  to  give 


Uq  =  —  exp{?'I'ir} 


T  S 


in  which  the  constant  k\  is  given  by 

l  _  ~kokij 

The  quantity  L  is  defined  by 


1  4 


(L*+2L(9*+k2)) 


}■ 


L  =  e2  +  kl  +  iM. 


(2.8) 


(2.9) 


(2.10) 
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A  few  comments  axe  in  order  before  completing  our  asymptotic  solutions  for  k  in 
the  large  wavenumber  limit.  First,  we  note  that  the  phase  function  9  determined  above  is 
constant  so  that  the  ‘fast’  dependence  of  the  perturbed  flow  could  equivalently  be  expressed 
in  terms  of  the  variable  r*  =  Re*r.  Thus  an  alternative  and  possibly  more  instructive 
method  to  solve  the  asymptotic  problem  for  k  would  be  to  use  a  multiple  scale  approach 

using  the  variables  Re  *r,  r,  Re~r r, _  Secondly,  we  notice  that  the  slow  dependence  of 

Uq  (and  hence  that  of  Vq,Wq  and  p^)  on  r  takes  the  usual  WKB  form  for  a  second  order 
differential  equation  in  that  it  is  proportional  to  {J  9(r)dr}z^' .  The  extra  exponential 
factor  is  introduced  because  the  differential  system  under  consideration  here  is  of  sixth 
order.  We  now  introduce  the  following  nomenclature: 

un  =  un  exp(i  Re1/2  9  r),  n  =  0,1,2,... 


with  similar  expressions  for  vn,  wn ,  and  p„. 

The  solution  u0  is  singular  in  the  limit  of  r  — »  0  so  there  must  be  an  inner  viscous 
boundary  layer  at  the  origin  in  order  to  smooth  out  the  singularity.  We  have  seen  that 
the  only  acceptable  roots  of  (2.7)  are  9\  and  —9\  since  the  complex  roots  will  lead  to 
exponential  growth  when  r  increases  or  decreases.  Thus  for  r  =  0(1),  the  appropriate 
forms  for  u0,u0  and  u'o  .determined  by  (2.6b)  and  (2.8)  are 


Lj  sin{Rei9jr  +  kjrj  _  NL\  cos{Rei9ir  +  k^r} 

— u0  _  — — - 1 - - - 

2rs  2 r» 

v  _  sin{Rei9ir  +  kir)  ^  N  cos{i?e i 9\r  +  k^r) 
r  s  r  * 

,  —  9\L\  cos{Ret9ir  4-  kir}  L\N9\  sin{i?ea0ir  +  jk.jr) 

—  KoWq  —  - T - (-  - r - 

2r\  2 ri 


(2-lla,  b,  c) 


Here  the  constant  L\  =  9\  4-  k%  +  iM,  while  N  is  a  constant  to  be  found  . 

We  can  see  from  (2.11)  that  we  cannot  choose  9\  and  N,  the  two  constants  at  our 
disposal,  in  order  to  satisfy  the  no-slip  condition  directly  at  r  =  1.  The  remedy  is  to  allow 
for  a  viscous  boundary  layer  at  r  =  1  in  which,  in  effect  ,  the  exponentially  decaying  WKB 
solution  can  be  used  to  satisfy  the  condition.  In  fact,  it  is  easiest  to  simply  use  a  boundary 
layer  variable 

£  =Re±{\ -r) 
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and  note  from  (2.11)  that  the  r  =  0(1)  core  solution  at  r  — ♦  1  can  be  expressed  as 
vo  — 1 >  sin  &i  (  [A'  sin(i?e  *  &x  4-  &i )  -  cos(J2e  ^B\  +  fcj )] 

i  -  ,  -  (212) 

+  costfif  [sin(i2e^i  4-  k\)  +  N  cos(iJe»  B\  +  fcj)] 
together  with  similar  expressions  for  uo,wo- 

In  order  to  match  with  the  above  limiting  form  of  the  core  solution,  we  must  expand 
( u,v,w,p )  in  the  £  =  0(1)  waill  layer  in  the  form 


(u,v  ,m,p)  =  (uo(?).i<o(f). 


(2.13) 


If  the  above  expansions  are  substituted  into  (2.1)  and  like  powers  of  Re  i  are  equated  we 


find  that 


^  ^2  —  ko  —  ^  “o  4-  2uo  =  —  Po 

{  ~  k°  ~  —  2u0  =  0 

{  ~  ~  ~  ~koPo 


-fio  +  k0w0  =  0. 


(2.14a) 


(2.146) 


(2.14c) 

(2.14d) 


If  we  eliminate  u0,u’o,Po,  we  find  that 


(jjp  -  *0  -  iM)  (j^S  -  *l)  -  ikl  «o  =  0 


(2.15) 


so  that  the  solution  for  Vo  that  satisfies  the  no-slip  condition  at  r  =  1  is 
€'o  =  a  [sin  B\(,  4-  A  cos  9i£  +  -  (A  +  B)e~,9i(] 


(2.16) 


Here  A,  B,  a  are  constants  to  be  determined  and  02, 63  axe  the  roots  of  (2.7)  that  have 
negative  imaginary  parts.  The  functions  uo,u>o,po  cam  then  be  found  from  (2.14)  and  the 
application  of  the  no- slip  condition  on  u0,u>o  shows  that 


B  =  — 


_ iL\6\ _ 

-L3B3  4-  L2B2  + 


A  B(L3  -  L2) 
Lx  -  L3 


(2.17a,  6) 
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while  the  matching  with  the  coreflow  requires  that 

a  =  N  sm(Re*8i  -f  ki)  -  cos(-Re*0i  4-  ki) 
Aot  =  sin(/?e*^i  4-  k\)  4-  N  cos (Re^di  4-  fci). 


Solving  these  equations  for  A  yields 

_  tan(.Re»fli  +  ki)  +  N 
N  tan(iJe  j  8\  4-  k\ )  —  1 


(2.17c) 


Since  A  is  already  known,  it  remains  for  us  to  determine  the  constant  N  and  then  (2.17c) 
will  constitute  an  eigenvalue  problem  for  ko  and  k\ .  This  constant  can  only  be  determined 
by  matching  the  coreflow  solution  with  the  solution  of  the  inner  boundary  layer  problem. 
2a  Inner  boundary-layer  problem 
Here,  we  define  r)  =  Re  *  r  and  write 


u  =  R€*u0(t /)  +  Re~*iii(r ))  +  .. . , 
v  =  Reivo(rj)  +  Re~*V!(Tf)  +  ..., 
w  =  Re*w0(rj)  +  Re~*wi(ri)  +  . . . , 
p  =  Re~*po(ri)  + .... 


(2.18) 


We  have  anticipated  in  (2.18)  that  the  disturbance  in  the  center  layer  is  Re*  larger  than 
in  the  core;  this  increase  in  size  is  necessary  because  of  the  dependence  of  the  core  solution 
on  r~  3 .  The  expansions  (2.18)  are  substituted  into  (2.1)  and  like  powers  of  Re~  *  are 
equated.  Surprisingly,  we  find  that  (u0,  Uo,  Wo,Po)  satisfy  the  full  equations  (2.1),  but  with 
r  replaced  by  rj,  k  by  fc0,  and  Re  =  1. 

Hence,  the  solution  cannot  be  completely  determined  analytically  in  the  core  so  (2.1) 
must  be  solved  numerically  with  Re  =  1.  However,  for  larger  values  of  tj  we  can  show  from 
(2.1)  with  Re  —  1  and  r  replaced  by  rj  that  a  large  rj  solution  of  (2.1)  is 


sin  #i77  a  cos  n  be~,$3,>  ce-1**17 

— r1  +  — +  — r—  +  — —■ 

r}i  rfi  TJ2  TJ* 


(2.19) 


together  with  similar  expressions  for  uq,xjuo  and  po-  Here  the  constants  a,b  and  c  can 
only  be  found  by  integrating  (2.1)  numerically.  This  was  done  by  first  calculating  three 
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independent  solutions  of  (2.1)  in  (0,Tjoo)  where  tj^  is  some  suitably  large  value  of  tj.  These 
three  solutions  are  then  multiplied  by  constants  such  that  (tt0,  v0,  tD0,po,  t^,  tBo)  determined 
in  this  way  is  continuous  at  with  the  same  function  determined  asymptotically  for  large 
r).  This  can  be  achieved  by  a  suitable  choice  of  these  three  constants  and  a,  6,  c,  in  (2.19). 
Clearly,  no  iteration  is  necessary  so  (2.1)  is  solved  numerically  only  three  times  in  this 
procedure.  The  matching  of  the  core  and  the  center  layer  solutions  then  requires  that 
N  =  a,  so  that  (2.17c)  yields 

Re*  +  A-i  =  nir  -f  tan-1  j  j  (2.20) 

which  determines  ko  and  k\.  In  fact,  since  Re  >>  1,  (2.20)  reduces  to 

i?e  *0i  =  nir,  h  =  tan"1  j  },  (2.21) 


where  n,  an  integer,  is  formally  O(Res).  We  can  take  to  be  defined  by  the  principal 
value  of  tan-1  [(AT  +  A)/(NA  —  1)]  since  adding  2 ir  to  ki  can  easily  be  shown  to  be  formally 
equivalent  (when  Re  >>  1)  to  a  different  choice  of  n  in  the  equation  for  9\. 

We  now  see  how  the  eigenvalues  k„  are  to  be  determined  for  Re  >>  1.  First,  we 
choose  a  large  integer  n  and  define 


*1  = 


nn 
Re  3 


Secondly,  we  solve  (2.7)  for  the  three  possible  values  of  k0  in  the  half  plane  R{fco}  >  0. 
For  each  of  these  values  of  ko,  we  solve  (2.1)  to  find  a  =  N,  then  fci  is  found  from  (2.9) 
and  (2.21). 

2b  A  comparison  with  exact  eigenvalues 

In  Figure  1  we  present  plots  of  the  imaginary  part  vs.  real  part  of  k  for  both  the 
present  approximations  and  the  exact  values  for  the  case  Re  =  1000,/  =  0.1.  The  two 
sets  of  data  are  seen  to  be  in  close  agreement. 

3  Some  observations  on  the  use  of  the  method 

1.  It  should  be  noted  that  the  method  is  formally  justified  when  the  core  variation  of 
the  disturbance  is  on  the  fast  Re  j  lengthscale.  In  effect,  this  means  that  the  eigenfunctions 
must  have  many  zeros  in  the  domain  of  interest.  Notwithstanding  this  observation,  in  the 
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Figure  1.  Plot  of  K[  vs.  KR  for  Re  =  1000,  f  =  0.1  (k  *  KR  +  i  Kf). 


spirit  of  asymptotics  we  cam  now  ignore  this  word  of  caution  and  apply  the  method  for 

the  case  n  =  1,2, 3, 4, 5, .  Though  it  is  formally  valid  for  the  higher  values  of  n,  we 

can  reasonably  expect  to  get  useful  answers  at  small  •values  too.  However,  since  the  formal 
requirement  for  the  validity  of  the  method  is  n  =  0(f?e»)  ,  the  method  could  well  be 
perverse  enough  to  give  better  guesses  for  the  lower  eigenvalues  at  the  lowest  values  of 
Reynolds  number  used. 

2.  The  order  in  which  the  eigenvalues  come  out  of  the  WKB  solution  is  related  to  the 
three  branches  we  have  discussed  elsewhere.  Essentially  the  eigenvalues  fall  into  groups 
of  three,  one  from  each  of  the  families.  No  sensible  comment  about  the  nature  of  the 
branch  relationships  for  the  first  few  eigenvalues  predicted  by  the  method  can  be  made 
since  inherently  the  method  is  being  used  out  of  its  range  of  validity;  all  that  one  can 
say  with  any  confidence  or  credibility  is  that  if  you  go  feu-  enough  down  the  series  the 
eigenvalues  predicted  by  the  WKB  method  will  fall  out  in  groups  of  threes,  one  from  each 
branch. 

3.  The  method  can  be  easily  modified  to  predict  eigenvalues  in  the  presence  of  an 
inner  rod.  There  are  really  just  two  cases  worth  considering.  If  the  rod  is  of  thickness 
Rt~ ,  the  above  formulation  holds,  but  the  inner  equations  now  have  to  be  solved  subject 
to  appropriate  boundary  conditions  at  a  scaled  value  of  r.  Much  more  realistically,  though, 
the  usual  situation  will  correspond  to  the  case  where  the  inner  boundary  is  at  a  finite  value 
of  r.  Suppose  then  that  there  is  such  a  boundary  at  r  =  q.  The  boundary  can  be  either  a 
free  surface  or  another  rigid  rod.  In  either  case  a  boundary  layer  of  thickness  Re~^  is  now 
set  up  at  the  inner  boundary  in  which  the  azimuthal  velocity  field  will  again  satisfy  (2.15) 
but  with  some  appropriate  variable  replacing  f.  A  solution  corresponding  to  (2.16),  but 
appropriate  to  the  inner  boundary  condition,  is  then  written  down.  Matching  with  the 
coreflow  then  gives  the  eigenrelation.  Now,  however,  there  is  no  inner  problem  to  be  solved 
numerically,  basically  because  r  is  not  small  enough  for  the  terms  like  £  to  be  comparable 
with  derivatives  with  respect  to  r  in  the  radial  direction. 
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Albuquerque,  NM  87110 

Aberdeen  Proving  Ground 

2  Commander,  CRDEC,  AMCCOM 
ATTN:  SMCCR-MU,  Mr.  W.  Dee 

SMCCR-RSP-A,  Mr.  M.  Miller 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes. 
Your  comments/answers  to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  BRL  Report  Number _ BRL-CR-630 _ Date  of  Report _ MAY  90 _ 

2.  Date  Report  Received _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest 

for  which  the  report  will  be  used.) _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source 
of  ideas,  etc.) _ 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars 
saved,  operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate 
changes  to  organization,  technical  content,  format,  etc.)  _ _ 


Name 


CURRENT  Organization 

ADDRESS  _ 

Address 


City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  New  or  Correct 
Address  in  Block  6  above  and  the  Old  or  Incorrect  address  below. 


Name 


OLD  Organization 

ADDRESS  _ 

Address 


City,  State,  Zip  Code 


(Remove  this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  mail.) 
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